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Abstract 

Four level quantum systems, known as quartits, and their relation to two- 
qubit systems are investigated group theoretically Following the spirit of 

2 ; 



Klein's lectures on the icosahedron and their relation to Hopf sphere fibra- 
tions, invariants of complex reflection groups occuring in the theory of qubits 

0\ 



o 



X 



and quartits are displayed. Then, real gates over octits leading to the Weyl 
group of E 8 and its invariants are derived. Even multilevel systems are of 



in 

interest in the context of solid state nuclear magnetic resonance. 
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1. Introduction 

A four-level system, also denoted a quartit, is a promising concept for the 
design of versatile two-qubit states and gates of quantum computation. As 



shown for instance in 



l|, states of nuclear spin | in a specific GaAs quantum 
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well device may be used for realizing the logical single and two-qubit gates by 
applying selective radio frequency pulses at the resonance frequency between 
two energy levels. Thus, quantum computing based on solid state nuclear 
magnetic resonance (NMR) with spins 3/2 (corresponding to quartits), spins 
7/2 (corresponding to 8-level systems, or octits) and higher order spins, is a 
strong motivation for our research. 

Multilevel systems with a number d of levels not a prime number (such as 
quartits, sextits and so on) are special by the mere fact that the number of 
mutually unbiased bases (MUBs) that can be constructed from the machinery 
of finite fields is strictly less than the value d + 1 corresponding to a complete 
set. It is well known that a complete set of MUBs may be obtained for a 
system of m qudits, if the qudit is a p-level system, i.e. d = p m and p a 
prime number 2] . It has been shown in previous works 3|, |4| that the free 
modules and the projective line over the finite ring Z2 x Z3 may be used to 
represent the geometry of commutation relations over a sextit system, leading 
to a maximum of three MUBs. In the context of MUBs, the sextit system 
is relevant because it corresponds to the smallest composite dimension for 
which a complete set cannot, in principle, be derived. But the sextit system 
cannot be ditinguished from the qubit/qutrit system [4] so that, from the 
point of view of quantum computation, it does not add new features. 

In contrast, in dimension four, a two-qubit system is not equivalent to a 
quartit system. A two-qubit system is obtained by taking tensor products 
of ordinary Pauli spin matrices a x , a y and a z , while the quartit system is 



generated by the two shift and clock operators (see Example 7 in |5| 
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Z = diag(l, u, u , ur) = cr z ® P, 



(1) 



where u; = exp(^ ZL ) and P = diag(l,i). 

The two-qubit Pauli group P 2 is generated by the two-fold tensor products 
of ordinary Pauli spin matrices. It is isomorphic to the small permutation 
gr oup [64, 266], the gr0U p of I1Um be r 266 in the sequence of s m al. gr oup 8 with 

cardinality 64. It may also be seen as a central product lj: V2 — Ef 2 * Z4 
since [64, 266] contain the extraspecial groups Ef 2 and the cyclic group Z 4 as 
normal subgroups, and Ef 2 H Z 4 coincides with the center Z 4 . 

The quartit group V qua rtit, generated by matrices X and Z, is isomorphic 
to the small group [64,18]. It may also be seen as a semidirect product: 
'P 'quartit = Z 4 x Z 4 ( or as a central product Z4 * Z 4 ). 

The reason why quartit actions relate to two-qubit actions is made ex- 



plicit in 



7|. In particular, the standard (entangling) cnot gate is easily 



inplemented by applying a suitable r.f. n pulse between two of the levels of 



the quartit system (see Sec. 3.3 of 



l|). This may be reformulated in the 



group language. The group C cno t = {X, Z, cnot) obtained from the quartit 



In this paper, symbols x, *, x and . denote the direct, central, semidirect and dot 
product of groups, respectively [6| 



group by adding to it the cnot generator reads C cno t — V% x (Z 2 x 54) (of 
order 3072). This means that two-qubit Pauli operators arise naturally in 
the normal subgroup part of C no t , as they do in the two-qubit Clifford group 
C 2 = V 2 x (Z 2 x S 6 ) (see (14) of [8]). The fourfold symmetry attached to 
the quartit system is visible in the factor group Ccnot/^2, that contains the 
symmetric group S4; similarly, the corresponding factor group for the two- 
qubit system contains the symmetric group Sq ( see also [9| about the sixfold 
symmetry of the two-qubit system). 

The connection of the quartit system to the two-qubit Clifford group is 
made more stringent by using (instead of the cnot gate) the entangling gate 
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that we introduced in several contexts 
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In some respect, gate S may be seen as a generalization of the single qubit 



Hadamard gate H = A= 



1 1 



\ 



V 



. The real gate S encodes the simulta- 



neous (entangled) eigenvectors of the triple {a x ® a x , a y ® a y , a z ® a z } [10] 
and lies in the automorphism group of the lattices Z 4 (see [2] in 11] ) and -D4 
(see [7] in [uj). 

The single qubit Clifford group is generated by H and P as C\ = (H, P) = 
V\ x Dq (with Dq the 12-element dihedral group) identifies to the rank 2 
complex reflection group Uq (number 9 in the Shephard-Todd sequence [i 
Similarly, it is straigthforward to check that the group Cg = (X, Z, S) 



V 2 x Sq, obtained by completing the quartit group with S (instead of cnot) 
is a subgroup of index two in C2, identifies to the rank 4 complex reflection 
group f/31 (number 31 in the Shephard-Todd sequence, see also p. 8, footnote 
gin U). 



2. Invariants involved in the theory of a single qubit 

There is a well known connection of invariant theory to Clifford groups 
and self-dual codes [12|, 113J. In the code context, the genus-n full weight 
enumerators of some self dual codes are polynomials in 2™ variables invariant 
under the Clifford group C n . Here, we are interested in a explicit connection 
of polynomial invariants of C\ and C$ to quantum states of the single qubit 
and quartits respectively In this goal, it is useful to recollect the celebrated 
Klein's lectures about the homogeneous invariants attached to the geometry 



of platonic solids 14J and give them an interpretation based on the Bloch 



sphere [15(|. Among Klein's invariants, fundamental invariants of C\ (of de- 
grees 8 and 24) are polynomials in the two amplitudes of a single qubit state, 
as shown in Sec. 12.21 The reasoning is generalized, in Sec. [31 to the fun- 
damental invariants of the complex reflection group C$ (of degrees 8, 12, 20 
and 24), which are polynomials in the four amplitudes of the quartit. 

2.1. Invariant ring of a finite matrix group 

Beforehand, let us recall some technicalities. An algorithm that calculates 
invariant rings of finite linear groups over an arbitrary field IK is described in 



Q 



171 ] and implemented in Magma 



6fl (see also [1J, U, 



18j). 



It may be applied to any finite matrix group G < GL n (K) of degree n 
over K. The group G acts linearly on the polynomial ring K[xi, . . . ,x n ] of 
the variables Xj q The invariant ring is defined as the set of polynomials left 
invariant under the action of G 



R = K[ Xl , . . .,x n f := {/ G K[ Xl , . . .,x n ]\a(f) = f, Va G G} . (2) 

Primary invariants may be constructed, i.e., homogeneous invariants /i, . . . , f n 
that are algebraically independant, such that the invariant ring is a finitely 
generated homogeneous module over A = K[/i, . . . , f n ]. 

Then, a (minimal) set of secondary invariants for R, with respect to these 
primary invariants, is a (minimal) generating set for R considered as a module 
over the algebra generated by the primary invariants. 

Given I a homogeneous ideal of the graded polynomial ring P = K[xi, . . . , x n ], 
then the quotient ring P/I is a graded vector space: P/I is a direct sum of 
the vector spaces Va for d — 0, 1, . . . where Vd is the K- vector space consisting 
of all homogeneous polynomials in P/I of weighted degree d. The Hilbert 
series of the graded vector space P/I is the generating fonction 



H P/I (t) = J2^MVd)t d , 

d=0 

and may be rewritten as a rational function of the variable t. 



2 In the sequel of this paper, variables Xi are interpreted as the complex amplitudes of 
a qudit and are denoted a, /3 (for a qubit) and a, /3, 7, S (for a quartit). 



If the character char(K) of K. does not divide the cardinality |G| of the 
group G (the non-modular case), the Hilbert series becomes the Molien series 
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2.2. Geometry of the invariants of the single qubit Clifford group C\ 

Let us see how the same results follow from the geometrical approach 
[141 ] . In the computational basis, a single qubit reads as the superposition 

|^) = «|0)+/3|1 ), a,(3eC, \a\ 2 +\f3\ 2 = l, (3) 

where the complex amplitudes a = a + ib, (3 = c + id satisfy the equations 
of the 3-sphere S 3 : a 2 + b 2 + c 2 + d 2 = 1. This state may also be represented 
(up to a global phase) in the Bloch sphere picture by using the map [l5| 

f = (<r x ) f = 2 Re(«/3), 
V = (oy)^ = 2 Im(a/3), 

■4, '— i ■ ! i ' i ■■ 

(4) 



C=(^L = I« 12 "' 2 



satisfying the equation of the 2-sphere 5* 2 : £ 2 + rf + ( 2 = 1. Specific points 
are the north pole |0 ), the south pole |1 ), and a set £ of four points of 
the equatorial plane located at the intersection of the Bloch sphere with the 
reference axes (±1,0,0) = ^(|0 ) ± |1 )), (0,±1,0) = ^=(i|0)±|l )). 

7 
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The Bloch sphere picture defines the Hopf fibration S 3 — » S 2 , in which 
the (great circle) fibre S 1 represents the global phase. 

The picture of the Bloch sphere may be supplemented by the picture of 
the Riemann sphere C U {00} 



f „,\ 



x 



\ V J 



Then, following Klein 
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\a\ 2 + \B\ 2 = l. 



(5) 



\ V J 



14j . points of the Riemann sphere minus the north 
pole (0, 0, 1) are mapped to the equatorial plane ( = 0, via a stereographic 
projection 



S(£,T),0 






(6) 



and the north pole is mapped to the fraction s(0, 0, 1) = i. 

In this representation, a point of the sphere S 2 can be seen as a equivalence 
class - (so that to fractions - and ^, such that pq' — p'q = 1, represent the 
same point on the sphere). For this reason, the Riemann sphere may also be 
identified to the complex projective line CP 1 . 

As a result, the fraction - of the Riemann sphere picture corresponds 
to the qubit a |0 ) + ft |1 ) of the Bloch sphere picture. The specific points 
above are the north pole g = |0 ), the south pole ° = |1 ), and points of the 
set £ that are s(±l,0,0) = ±j = ^(|0 ) ± |1 }) and s(0,±l,0) = ±{ = 
^(|0>±i|l». 

Having parametrized the 2-sphere as CP 1 , one can capture the symmetry 



of a set of n points S = <— ,■■■, — >, living on it, in a polynomial of the 
form pj, [l6| 

&- V -±)-&- V -), (7) 

P qi p q n 



that has roots exactly at the points of S. This can be rewritten as a homo- 
geneous polynomial of degree n 

(q 1 a-pij3)---(q n a-p n j3). (8) 

The pairs (a, /3) which cancel the polynomial (jSJ) are the amplitudes of the 
quantum states in ([3]). For instance, the set of six points listed above may be 
seen as the vertices of a octahedron and their symmetry may be expressed 
as the polynomial 

T := a/3(a -0)(a + f3)(a + i/3)(a - t/3) = a/3(/3 4 - a 4 ), (9) 

which is invariant under the action of the octahedral group. 

The tetrahedral group 

Refering to the specific points displayed at the previous section, the faces 
of the octahedron are centered at the vertices of a cube of coordinates 



±Z = ±r l = ±C = ^=, (10) 



and one can pick up the vertices of a tetrahedron by selecting, among the 
eight possible combinations of sign, those for which the product £r)( is posi- 
tive 

p 1 + 2 1 — i —1 + i —1 — i 



q V3-l' V3 + 1' V3 + 1' V3-l' 

A straightforward calculation leads to the invariant polynomials 



(11) 



which correspond to 



a 4 ± 2iV3a 2 (3 2 + /3 4 , (12) 



a tetrahedron (upper sign) or to a countertetrahedron 



(lower sign) (see 14j, p 54). One can check that the tetrahedral group 



U\ = SL(2,3) (number 4 in the Shephard-Todd sequence) 



1/ 1 — i 1 — i 11/ 1 — i i — 1 

- ' -v^ -1 1 + VZ \ 1 + V3 l + \/3 



possesses an invariant ring spanned by the two primary invariants ([H]) and 
( fl2|) . corresponding to the Molien Series 



MS(Kt) = 7 TT7 ^r- (14) 

V ; (l-t 4 )(l-t 6 ) V ; 



The tetrahedral system is involved in the construction of minimal four- 
state quantum tomography. The quantum states correspond to the four 



vertices of the tetrahedron and form a SIC POVM (see eq (2.6) in [20(] and 



10 



Sec. B in 21]). They may be taken as 



\if>) = a\Q)+p\l), a x \ij>), cT y \^), (T g \il>) 



with a = -Wl + l/v^) 172 and /3 = -iexp(z7r/4)(l - l/v^) 172 , 



(15) 



in agreement with the upper equation ( 1121 . 



T7ie octahedral group 

The octahedral group is not of the reflection type but is isomorphic to C/4 



22|. It possesses the invariant polynomial ([9]) for the vertices and another 
invariant polynomial related to the 8 faces of the octahedron centered at the 
points given in ffTUj) . It is obtained by multiplying together the invariants 
(EE2D so that 

W := a 8 + Ua 4 (3 4 + /3 8 . (16) 

The octahedral group O is generated as the derived subgroup C[ of the single 
qubit Clifford group 



0-{irr t l\ l % j }. (17) 

- l 1+i l + i 



The Molien series is 

1 - t 4 + t 8 



MS(O) 



(l-t6)(l- t B) 

and the invariant ring is spanned by (Q and (1X6]) . 
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The group Us and the Clifford group C\ 

The reflection group Us — Z±. S^ is a subgroup of index two in the Clifford 
group C\ 




(19) 



Its Molien series is 



MS(J7 8 



(20) 



(l-£ 8 )(l-£ 12 )' 

and the invariant ring is spanned by the invariant W and an invariant of 
degree 12 

K := « 12 - 33 a 8 /? 4 - 33 /3V + /3 12 . (21) 

There exists the following relationship (eq (53) in 14 1) |f| 



108 T 4 - W J + k z = 0. 



(22) 



Finally, the Clifford group C\ is the reflection group U$. The Molien series 
1/(1 — t 8 )(l — t 24 ), is spanned by the three invariants T, W and k. 

In 22|, the invariant ring C[a, /5] Cl of C\ is found to be isomorphic to the 
polynomial ring C[i?4, A 12 ] generated by the Eisenstein series E4 of weight 
4 and the cusp form A 12 of weight 12, but no reference to Klein's work is 



3 Primary invariants T and W, the secondary invariant k and relation (|22[) follow 
from Magma with the following code "R:=InvariantRing(U8); Primarylnvariants(R); Sec- 
ondarylnvariants(R); Algebra(R); Relations(R);" 



12 



pointed out. The invariant ring of C\ is spanned by the weight enumerator 
W of the Hamming code e% and the weight enumerator Q of the Golay code 



G 24 (12 1 



g : = a 24 + 759 a 16 /? 8 + 2576 a 12 /? 12 + 759 a 8 ft 16 + ft 24 , (23) 



where T 4 



42 • 



3. Invariants involved in the theory of a single quartit 

Let us define a quartit state as 

\ij)=a\0)+(3\l )+ 7 |2) + <J|3), (24) 

with a, ft, 7, 5 G C and the sum of squared amplitudes equal to unity. 

Formally, one can map a quartit state to a two-qubit state as |0 ) = |00 ), 

|1 ) = |01 ), |2 ) = 1 10 ) and |3 ) = |11 ), so that 

\iP ) = a 1 00 ) +0 |01 ) +7|10 } + S |11 ) , (25) 

with a, ft, 7, S E C and the sum of squared amplitudes equal to unity. This 
key identification is carried out in [l]. Although a quartit is not, strictly 
speaking, equivalent to two qubits, both systems are close to each other as 
soon as the suitable gate action is performed (as shown in the introduction). 
Experimentally, it is much convenient to have at one's disposal a four-level 

13 



system (such as the NMR spin | state) than a two-qubit system. In the 
present context of the second Hopf fibration, it is true that a quartit is for- 
mally similar to a two-qubit system, and the resulting polynomial invariants 
also are the same. 

The quantum states are now the inhabitants of the 7-sphere S 7 of unit 
radius. The latter can be mapped to the 4-sphere S A by using , instead of 
complex numbers, the quaternions 

Qi = a + (3j, Q 2 = 7 + Sj, with Q u Q 2 e H, (26) 

where EI denotes the (non-commutative) field of quaternions. The multipli- 
cation rules are i 2 = j 2 = k 2 = — 1, ij = k, jk = i, ki = j, and ji = —k, 
kj = —i and ik = —j. The conjugate of a quaternion Q\ is Qi = ot — (3j and 

the norm is |Qi| 2 = QiQi- 

s 3 
In the second Hopf fibration S — > S , the coordinates of a point of the 

target space are (£, r],u, v,Q — { a i)^ {i — 1-5) where 



. 1 \ / i,j,k \ / 1 

CTl = <T X = I , cr 2 .3,4 = , &5 = O z = | , (27) 

10/ V -(*.i.*) ° / V ° - 1 



are quaternion Pauli matrices 15| . Explicitely 



^ = 2Re(a7 + /35), n = 2 Im(a7 + 05), 
u = 1 Re(a5- j 07), u = 2 Im(a5 - ^7), 

C = IQ?|-|Q^I = « 2 -/3 2 + 7 2 -5 2 . 

14 



(28) 



The notable point is the sensitivity of the second Hopf map to entanglement. 
Separable states satisfy the condition aS = f3 r y and therefore are mapped onto 
a subset u = v — of pure complex numbers in H. Conversely, the maximally 
entangled states correspond to a subset of pure quaternions £ = r\ = ( = 
inH. 

It is known that the projective Hilbert space of two non-entangled qubits 
is the product of two 2-dimensional spheres S 2 , each sphere being the Bloch 
sphere attached to a specific qubit. In terms of the second Hopf fibration, one 
sees that the space S A maps to a unit sphere S 2 (such that u — v — 0) which 
is the Bloch sphere for the first qubit. The second Bloch sphere is recovered 
from the fibre. The projective Hilbert space for maximally entangled states 
is known to be S* 3 /^, i.e. a 3 sphere such that two opposite points are 
identified. This has a counterpart in the Bloch sphere picture in the fact 
that opposite points on S 3 corresponds to the same maximally entangled 



state [15] 



The Riemann sphere CUoo may be generalized to the one-point com- 
pactification HUoo of the quaternion numbers, which may be identified to 
the projective line HP 1 = S A over the quaternions. To proceed further, one 
can map points of HP 1 to the plane ( = 0, via a 4-dimensional stereographic 



15 



projection 



*(£, V , u, v, = ^ + %V +^ + - with S (0, 0, 0, 0, 1) = i (29) 



but the generalization of Klein's approach is not straigthforward since a 
homogeneous polynomial corresponding to (J7|) should display four variables 
a, (3, 7 and S instead of two. Fortunately, the invariant theory of finite linear 
groups summarized in Sec. 12.11 may be applied and results made explicit by 
using Magma 6]. 

The invariant ring of complex reflection group C$ = U$x attached to a 
quartit is spanned by invariants of degrees 8, 12, 20 and 24 in the four 
variables x% = a, x 2 = (3, x 3 = 7 and x 4 = 5. The smallest degree invariant 
is 

inv 8 := S 8 + 14 £4,4 + 168 £ 2)2 , 2i2 , (30) 

in the notations of [12|,ll3j, i.e. S 8 = ^=1 x h ^4,4 = Ylj>i x \ x ) an d S 2j2)2j2 = 
Invariant inv§ also represents the genus-two complete weight enumerator 



of the code e 8 ® F 4 12|, and indeed generalizes Klein's invariant ( 1T6|) . 



With the same type of notations, higher order invariants are as follows 

inv 12 := E12 - 33 S 8j4 + 330 £ 4A4 + 792 E 2 , 2 ,2,6, 

inv 20 := S 20 - 19 E 16)4 - 494 E 12]8 + 380 E 12 , 4 , 4 + 7296 E 10 , 6)2j2 

+1710 E 8jM + 133380 E 8)4 , 4) 4 + 102144 E 6)6j6)2 , 

16 



inv 24 := £ 24 + 759 E 16j8 + 2576 E 12jl2 + 212520 E 12A4 ,4 + 340032 S 10 , 6 , 6 , 2 
+22770 S 8i8>8 + 1275120 £ 8jM)4 + 4080384 S 6i6j6i6 . 

(31) 

Note that inv 12 and inv 24 generalize k in ( 121 j) and Q in (j23p . and that inv 24 
is the genus-2 Hamming weight enumerator of the Golay code G 24 . 

4. On the Weyl group of E s generated from octits 

A relevant example of an eight-level (spin |) system, here also denoted a 



octit, is the 133 Cs in an anisotropic environment 23] . Real gates over a octit 
are especially interesting because they allow to generate the largest finite 
complex reflection group Use — W(E$), of order 696729600, where W(E 8 ) 
denotes the Weyl group of Lie algebra E 8 (see also [lOj, [llj for a similar 
approach of Es). One obtains 

^ 36 = <^8,/®/®^,5 3 ), (32) 

where I is the 2x2 identity matrix, X 8 is the 8x8 shift matrix and S3 
is a generalization of matrix S, that encodes the eigenstates of the triple of 
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observables o- z ®{a x ®a x ,ay®a yi a z ®a z } 



S* = 
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(33) 



The Molien series for Use reads 



n 



m=2,8, 12, 14,18,20,24,30 



l-t r ' 



(34) 



leading to invariants of degrees given by the indices in the product. The 
invariants of smallest degree are found as 



12 — S2, Is — Sg + 56 ^4,2,2 — 42 S4 ; 4 — 168 2^2,2, 2, 2- 



(35) 



Higher order invariants could not be obtained with this method on MAGMA, 
due to a lack of memory, even on a 96 MB segment of the cluster at our 



university. For recent work on this subject, see 271 ] . 
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5. Conclusion 

Complex reflection groups of the Shephard-Todd sequence happen to be 
the natural players in the context of the multilevel approach of quantum 
computation. Although their construction is well documented in mathemat- 
ics and the coding theory of self dual codes, due to the identification of the 
invariants to the complete weight enumerators, their introduction for ap- 
proaching the design of qubits, quartits and octits seems to be novel. The 
largest complex reflection group W(E$) naturally appears from the relevant 
real gates acting over a octit. It is fascinating that fingerprints of the radii 
of the Gosset dries over E% already appeared in an experimental solid state 



system (Ising chain) 



24 



25| . The Gosset polytope of the E 8 lattice may be 



embedded in the Hopf sphere S 7 of a two-qubit system, as shown in 



26|. 



The organic relationship of even Euclidean lattices to entanglement is also 
111]- 



explored in 



References 

[1] Hirayama Y, Miranovwicz A, Ota T, Yusa G, Muraki K, Ozdemir S K 
and Imoto N 2006 Nanometer-scale nuclear-spin device for quantum 
information processing. J. Phys.: Condens. Matter IS S885. 

[2] Planat M, Rosu H and Perrine S 2006 A survey of finite algebraic gro- 
metric structures underlying mutually unbiased quantum measurements. 
Founf. Phys. 36 1662. 



19 



[3] Planat M, Baboin A C and Saniga M 2008 Multi-line geometry of 
qubit/qutrit and higher order Pauli operators Int. J. Theor. Phys. 47 
1127. 

[4] Havlicek H and Saniga M 2008 Projective ring line of an arbitrary single 
qudit J. Phys. A: Math. Theor. 41 015302. 

[5] Kibler M 2009 An angular momentum approach to quadratic Fourier 
transform, Hadamard matrices, Gauss sums, mutually unbiased bases, 
the unitary group and the Pauli group J. Phys. A: Math. Theor. 42 
353001. 

[6] Bosma W, Cannon J and Playout C 1997 The Magma algebra system 
I. The user language J. Symb. Comp. 24 235. 

[7] Kessel A R and Ermakov V L 1999 Multiqubit spin JEPT Letters 70 
61. 

[8] Planat M and Kibler M 2010 Unitary reflection groups for quantum fault 
tolerance J. Comp. Theor. Nanoscience 7 1. 

[9] Planat M and Saniga M 2008 On the Pauli graphs of iV-qudits Quant. 
Inf. Comp. 8 127. 

[10] Planat M 2010 Three-qubit entangled embeddings of CPT and Dirac 
groups within E8 Weyl group Int. J. Theor. Phys. 49 1044. 



20 



[11] Planat M 2010 Entangling gates in even Euclidean lattices such as the 
Leech lattice. Preprint 1002.4287 [quant-ph]. 

[12] Nebe G, Rains E M and Sloane N J A 2001 The invariants of the Clifford 
groups Designs, Codes and Cryptography 24 99. 

[13] Nebe G, Rains E M and Sloane N J A 2006 Self-dual codes and invariant 
theory (Springer: Berlin). 

[14] Klein F Lectures on the icosahedron and the solution of equations of the 
fifth degree (Dover: New- York). 

[15] Mosseri R and Dandoloff R 2001 Geometry of entangled states, Bloch 
spheres and Hopf fibrations J. Phys. A: Math. Gen. 34 10243. 

[16] Weeks J 2006 Exact polynomial eigenmodes for homogeneous spherical 
3-manifolds Class. Quant. Grav. 23 6971. 

[17] Kemper G 1996 Calculating Invariant Rings of Finite Groups over Ar- 
bitrary Fields J. Symb. Comp. 27 171. 

[18] King S 2007 Minimal generating sets of non-modular invariant rings of 
finite groups. Preprint 0703035 [math- AC]. 

[19] Kane R 2001 Reflection groups and invariant theory (Berlin: Springer). 

[20] Rehacek J, Englert B-G and Kaszlikowski D 2004 Minimal qubit tomog- 
raphy Phys. Rev. A 70 052321. 



21 



[21] Durt T 2008 Wigner tomography of two-qubit states and quantum cryp- 
tography Phys. Rev. A 78 042338. 

[22] Bannai E 1999 Modular invariance property of association schemes, type 
II codes over finite rings and finite abelian groups and reminiscences of 
Franois Jaeger (a survey) Ann. Inst. Fourier, Grenoble 49 763. 

[23] Khitrin A K and Fung B M 2001 NMR simulation of an eight-state 
quantum system Phys. Rev. A 64 032306. 

[24] Coldea R, Tennant D, Wheeler E M, Wawzrzynska E, Prabhakaran D, 
Telling M, Habicht K, Smeibidl P and Kiefer K 2010 Quantum criticality 
in an Ising chain: experimental evidence for emergent E 8 symmetry 
Science 327 177. 

[25] Kostant B 2010 Experimental evidence for the occurence of E 8 in nature 
and the radii of the Gosset circles Select. Math., New Ser. 16 419. 

[26] Mosseri R and Ribeiro P 2007 Entanglement and Hilbert space geometry 
for systems with a few qubits Math. Struct, in Comp. Science 17 1117. 

[27] Talamini V 2010 Flat bases of invariant polynomials and P-matrices of 
E 7 and E 8 . Preprint 1003.1095 [math-ph]. 



22 



